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ABSTRACT — The  goai  of  this  study  was  to  determine  the  effectiveness  of  using  traditionai  and  new  vitai  signs  (heart  rate 
variabiiity  and  compiexity  [FIRV,  HRC])  for  predicting  mortaiity  and  the  need  for  iife-saving  interventions  (LSIs)  in  prehospitai 
trauma  patients.  Our  hypothesis  was  that  statisticai  regression  modeis  using  traditional  and  new  vital  signs  would  be 
superior  in  predictive  performance  over  models  using  standard  vital  signs  alone.  This  study  involved  108  prehospital  trauma 
patients  transported  from  the  point  of  injury  via  helicopter.  Heart  rate  variability  and  HRC  were  calculated  using  criterion 
standard  R-R  interval  sequences  manually  verified  from  the  patients’  electrocardiograms.  Means  and  standard  deviations 
for  vital  signs,  HRV,  HRC,  and  Glasgow  coma  scale  (GCS)  scores  were  obtained  for  nonsurvivors  versus  survivors  and  LSI 
versus  non-LSI  patient  groups  and  then  compared  using  Wilcoxon  statistical  tests.  Receiver-operating  characteristic  curves 
were  also  obtained  to  compare  different  regression  models  for  predicting  mortality  and  the  need  for  LSIs.  Seventeen 
patients  (16%)  died.  Eighty-two  patients  (76%)  received  a  total  of  142  LSIs.  Receiver-operating  characteristic  curves 
demonstrated  better  prediction  of  mortality  and  LSI  needs  using  heart  rate  and  HRC  (area  under  the  curve  [AUC];  AUCs, 
0.86  and  0.86)  than  using  heart  rate  alone  (AUCs,  0.79  and  0.57).  Likewise,  receiver-operating  characteristic  curves 
demonstrated  better  prediction  using  total  GCS  score  and  HRC  (AUCs,  0.82  and  0.97)  than  using  total  GCS  score  (AUCs, 
0.81  and  0.91).  Similar  results  were  obtained  for  heart  rate  and  HRV  (AUCs,  0.86  and  0.73).  The  major  implication  of  this 
study  was  that  traditional  and  new  vital  signs  (HRV  and  HRC)  should  be  used  simultaneously  to  improve  prediction  of 
mortality  and  the  need  for  LSIs  in  prehospital  trauma  patients  during  all  echelons  of  trauma  care.  Improvements  in  the  timely 
use  and  diagnostic  accuracy  of  transportable  vital  signs  monitors  will  require  use  of  traditional  and  new  vital  signs  from  the 
trauma  patient  cohort. 
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INTRODUCTION 

Initiation  of  early  and  effective  life-saving  interventions  (LSIs) 
is  an  important  aspect  of  trauma  medicine,  especially  in  the 
prehospital  and  battlefield  settings.  In  these  environments,  po¬ 
tential  delays  in  administering  LSIs  could  lead  to  increased 
morbidity  and  mortality  after  trauma  (1,2).  Traditionally,  medics 
have  relied  on  standard  vital  signs  such  as  blood  pressure  and 
heart  rate  (HR)  to  monitor  trauma  patients  for  abrupt  changes 
that  indicate  the  need  for  LSIs.  Unfortunately,  there  are  limi¬ 
tations  inherent  in  monitoring  traditional  vital  signs  that  could 
possibly  prevent  the  medic’s  ability  to  identify  patient  desta¬ 
bilization  until  late  and  irreversible  changes  in  state  take  place 
(1,  3-7).  Because  monitoring  of  standard  vital  signs  has  be¬ 
come  routine  throughout  the  spectrum  of  critical  care,  earlier 
administration  of  LSIs  implies  a  need  for  improvements  in 
the  timely  use  and  diagnostic  accuracy  of  vital  signs  moni¬ 
tors  (1,  2,  8). 
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One  possible  solution  for  improving  LSI  performance  via 
monitoring  may  be  the  use  of  new  vital  signs  such  as  advanced 
indices  derived  from  the  electrocardiogram  (ECG),  namely, 
HR  variability  (HRV)  and  complexity  (HRC),  which  are  ca¬ 
pable  of  detecting  subtle  changes  in  HR  behavior  driven  by  the 
autonomic  nervous  system,  which  could  precede  changes  in 
standard  vital  signs  (9-13).  These  indices  could  be  integrated 
into  existing  vital  signs  monitors  without  adding  new  sensors 
and/or  boxes  to  the  medic’s  kit.  Previous  studies  have  shown 
the  predictive  power  of  HRV  and  HRC  for  risk  stratification, 
diagnosis,  and  continuous  monitoring  of  traumatically  injured 
patients  (14).  Benefits  of  using  HRV  and  HRC  include  their 
ease  for  rapid  calculation,  noninvasive  nature,  and  robustness 
against  certain  kinds  of  noise  (15,  16). 

However,  as  standard  vital  signs  alone  may  not  be  reliable 
in  many  situations,  HRV  and  HRC  have  several  significant 
drawbacks,  including  limited  use  in  the  presence  of  high  arti¬ 
ficial  noise  levels  and  ectopic  beats  in  the  captured  waveforms 
(14—17).  Recent  work  suggests  that  new  vital  signs  such  as  HRV 
and  HRC  should  be  used  in  conjunction  with  traditional  vital 
signs  to  achieve  more  accurate  diagnostic  capabilities  (2).  In  light 
of  this  work,  the  objectives  of  this  study  were  1)  to  confirm  that 
HRV  and  HRC  can  discriminate  between  nonsurvivors  and  sur¬ 
vivors  of  trauma  and  2)  between  trauma  patients  who  received  one  or 
more  LSIs  and  those  who  received  none;  and  3)  to  compare  different 
multivariate  logistic  regression  models  for  predicting  mortahty  and  4) 
for  predicting  the  need  for  LSIs  in  trauma  patients.  Our  hypothesis 
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was  that  models  incorporating  traditional  and  new  vital  signs  to 
predict  mortality  and  the  need  for  LSIs  would  be  able  to 
outperform  models  incorporating  only  traditional  vital  signs. 

MATERIALS  AND  METHODS 
Subjects  and  protocol 

This  Study  was  conducted  under  a  protocol  reviewed  and  approved  by 
the  US  Army  Medical  Research  and  Materiel  Command  Institutional  Review 
Board  and  in  accordance  with  the  approved  protocol.  After  approval,  108  pre¬ 
hospital  patient  records  were  selected  from  the  US  Army  Institute  of  Surgical 
Research  Trauma  Vitals  (TV)  database  based  on  availability  of:  1  -  vital  signs 
data;  2  -  ECG  waveforms;  and  3  -manual  verification  of  all  R-to-R  interval  (RRI) 
sequences.  Because  all  data  were  analyzed  post  hoc,  the  study  was  considered 
minimal  risk,  and  infonned  patient  consent  was  waived. 

The  cohort  used  in  this  study  consisted  of  prehospital  trauma  patients  with 
blunt  and  penetrating  injuries  transported  from  the  point  of  injury  by  helicopter 
service  to  a  level  I  trauma  center  in  Houston,  Texas  or  San  Antonio,  Texas.  All 
patients  were  monitored  during  transport  using  a  Welch  Allyn  PIC  50  (Welch 
Allyn,  Skaneateles  Falls,  NY)  monitor.  Data  from  a  personal  digital  assistant 
and  flash  cards  attached  to  the  monitor  were  extracted  by  research  person¬ 
nel  and  uploaded  to  the  TV  database  for  analysis.  All  nonelectronic  data  were 
manually  recorded  on  the  run  sheet  from  the  monitor’s  screen  by  Emergency 
Medical  Services  medics,  then  collected  on  a  standardized  form  and  entered 
into  the  TV  database.  These  included  demographic  data,  physical  examination 
results,  Glasgow  Coma  Scale  (GCS)  scores,  and  interventions  performed  on 
the  patients  in  the  field.  Life-saving  inteiwentions  consisted  of  endotracheal 
intubations,  transfusions,  tube  thoracostomies,  cardiopulmonary  resuscitations, 
needle  decompressions,  angioembolizations,  cricothyrotomies,  thoracotomies, 
and  cardioversions. 

Importantly,  numeric  data  were  stored  at  a  rate  of  1  measurement  every 
3  min,  coinciding  with  the  patient’s  noninvasive  blood  pressure  measurements, 
whereas  wavefoiTn  data  were  stored  at  a  rate  of  375  Hz.  Of  the  108  wavefoiTns, 
17  ECGs  (16%)  belonged  to  nonsurvivors,  whereas  91  ECGs  belonged  to 
survivors  of  trauma.  In  addition,  82  ECGs  (76%)  belonged  to  patients  who 
received  at  least  one  LSI,  whereas  the  remaining  26  ECGs  (24%)  belonged  to 
patients  who  received  none.  Lengths  of  patients’  ECGs  varied  from  approxi¬ 
mately  15  to  20  min.  Heart  rate,  systolic  blood  pressure  (SBP),  diastolic  blood 
pressure  (DBP),  mean  arterial  pressure  (MAP),  respiratory  rate  (RR),  blood 
oxygenation  (Sp02),  shock  index  (SI  =  HR/SBP),  and  pulse  pressure  (PP  = 
SBP  -  DBP)  were  also  captured  during  trauma  care  for  patient  assessment. 

Heart  rate  variability  and  complexity 

For  this  study,  HRV  and  HRC  were  derived  using  true  RRI  sequences 
extracted  from  the  patients’  ECG  waveforms.  The  criterion  standard  for  obtaining 
true  RRI  sequences  was  manual  verification  of  R  waves,  which  was  accom¬ 
plished  by  importing  the  ECG  waveform  data  into  WinCPRS  software  (Abso¬ 
lute  Aliens  Oy,  Turku,  Finland),  visually  analyzing  the  data,  and  marking  times 
and  positions  of  all  R  waves.  An  HRV  ratio  was  computed  using  the  Poincare 
variability  ratio,  SD1/SD2,  which  describes  the  short-teiTn  and  long-term  RRI 
variability  of  a  sequence.  This  ratio  was  selected  owing  to  its  suitability  for 
analysis  of  shorter  time  series  and  frequent  citations  in  literature.  SDl  and  SD2 
were  obtained  from  the  following  equations: 

SDC  =  var  ,  1 

\\/2  \J2  /  2 

and 

5Z)2^  -  2^DRR^  -  ,  2 

2 

V/  =  1,  ...,  V,  where  var(«)  denotes  the  variance  function,  SDRR  denotes 
the  standard  deviation  of  RRIs,  SDSD  denotes  the  standard  deviation  of  suc¬ 
cessive  RRI  differences,  and  N  denotes  the  number  of  RRIs  in  the  sequence 
window  (18).  An  HRC  value  was  computed  via  the  method  of  sample  entropy, 
SampEn  (w,  r,  V),  which  describes  the  conditional  probability  that  2  epochs 
similar  for  m  RRIs  remain  similar  at  the  next  RRI,  i.e.,  differ  by  no  more 
than  some  tolerance  r  (in  milliseconds)  (19).  Parametric  values  (V  =  200, 
m  =  2,  r  =  6)  were  established  from  previous  work  (2,  12,  16,  17).  In  addition, 
the  following  equations  were  used  to  calculate  SampEn: 

SampEn{m,/-,  V)  —  -  In  [A/B)  ,  3 

B  =  [(A'-m-l)/2]l"7'B;.-(m)  ,  4 

A  =  [(iV-m-l)/2]X"-X(m)  .  5 

Here,  xj^i)  denotes  a  segment  of  m  consecutive  RRIs  starting  at  index  /  and 
running  from  i=  1 ,  denotes  the  number  of  epochs  x^ij)  within 


r  of  A'„,(/),  for  i  j,  multiplied  by  (V  -  w  -  1)  \  and  denotes  the  number 

of  epochs  Xm+iif)  within  r  of  x„j+j{i),  for  i  ^  j,  multiplied  by  {N  -m-  \)  ^  (19). 

Statistical  analyses 

Like  the  statistical  analyses  described  in  (2),  all  data  sets  were  analyzed 
using  Wilcoxon  tests  for  nonparametric  distributions.  Data  were  expressed 
as  means  ±  SD.  The  power  of  demographics,  vital  signs,  HRV,  HRC,  and  GCS 
scores  to  identify  nonsurvivors  and  whether  LSIs  were  performed  was  esti¬ 
mated  using  multivariate  logistic  regression  modeling.  In  other  words,  initial 
multivariate  logistic  regression  analyses  were  done  for  all  subjects  with  inde¬ 
pendent  variables  of  age,  height,  race,  and  weight  and  with  dependent  variables 
of  HR,  SBP,  DBP,  MAP,  RR,  and  SI.  These  analyses  excluded  HRV  and  HRC 
values.  Factors  that  were  not  significant  {P  >  0.05)  were  removed  from  the 
model  via  backward  elimination.  A  second  set  of  analyses  was  done  for  de¬ 
pendent  variables  of  HR,  SBP,  DBP,  MAP,  RR,  and  SI,  including  HRV  and 
HRC  for  perfoiTnance  comparisons  with  the  initial  set.  In  addition,  a  third 
set  and  a  fourth  set  of  analyses  were  performed  for  all  subjects  to  include 
GCS  scores,  with  and  without  HRV  and  HRC  as  dependent  variables,  re¬ 
spectively.  The  ability  of  statistical  models  to  predict  mortality  and  the  need 
for  at  least  one  LSI  was  assessed  using  sensitivity  and  specificity  scores  as  well 
as  receiver-operating  characteristic  (ROC)  curves.  Statistical  analyses  were 
performed  using  IMP  version  9.0.0  (SAS  Institute,  Cary,  NC)  and  the  R  Lan¬ 
guage  (http  ://w  w  w  .r-projec  t.  org/) . 

RESULTS 

Patients’  demographics  for  108  patients  are  tabulated  in 
Table  1,  whereas  interventions  performed  on  these  patients  and 
classified  as  life  saving  by  a  multidisciplinary  team  of  trauma 


Table  1.  Demographics  of  patients  (n  =  108) 


Characteristics 

Mean  ±  SD  or  n  (%) 

Age,  yrs 

37  ±  14 

Sex 

Male 

82  (76) 

Female 

25  (23) 

Not  recorded 

1  (1) 

Race 

White 

44  (41) 

Black 

6(6) 

Hispanic 

24  (22) 

Asian/Pacific 

3(3) 

Not  recorded 

31  (28) 

Mechanism  of  injury 

Blunt 

93  (86) 

Penetrating 

13  (12) 

Not  recorded 

2(2) 

Receiving  LSI(s) 

82  (76) 

Died 

17(16) 

Injury  severity  score 

17  ±  11 

Total  GCS  score 

9  +  5 

Field  HR,*  beats  per  minute 

106  ±  27 

Field  SBP,*  mm  Hg 

116  +  23 

Field  DBP,*  mm  Hg 

79  +  16 

ED  HR,*  beats  per  minute 

98  +  25 

ED  SBP,*  mm  Hg 

114  +  27 

ED  DBP,*  mm  Hg 

75  +  21 

‘Denotes  entry  values  taken  from  the  run  sheet. 

SD,  standard  deviation;  WVSM,  wireless  vital  signs  monitor;  ED,  emer¬ 
gency  department. 
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Table  2.  Life-saving  interventions 


Life-saving  interventions 

#  n 

%  n/142 

Prehospital 

87 

61 

Emergency  department 

55 

39 

Angioembolization 

1 

1 

Blood 

30 

21 

Cardiopulmonary  resuscitation 

9 

6 

Chest  tube 

18 

13 

Intubation 

73 

51 

Needle  decompression 

8 

6 

Surgical  cricothyrotomy 

2 

1 

Thoracotomy 

1 

1 

Total 

142 

100 

experts  are  listed  in  Table  2.  Of  these  108  patients,  17  (16%) 
died,  of  which  16  (15%)  required  an  LSI.  Of  the  108  patients, 
26  (24%)  did  not  require  an  LSI.  The  other  82  (76%)  patients 
received  a  total  of  142  LSIs.  Eighty-seven  (61%)  of  the  LSIs 
were  performed  prehospital  and  55  (39%)  in  the  emergency 
room.  Interventions  consisted  of  the  following:  73  endotracheal 
intubations,  30  transfusions,  18  tube  thoracostomies,  9  cardio¬ 
pulmonary  resuscitations,  8  needle  decompressions,  1  angioem- 
bolization,  2  cricothyrotomies,  and  1  thoracotomy.  Importantly, 
the  demographics  of  the  chosen  population  included  HRs 
ranging  from  20  to  156  beats  per  minute,  SBPs  ranging  from 
80  to  170  mm  Hg,  DBFs  ranging  from  44  to  1 10  mm  Hg,  and 
various  types  of  injuries  and  LSIs.  This  cohort  provided  the 
ECG  morphology  for  HRV  and  HRC  calculations. 

Means  of  HRV  and  HRC  statistics,  SDs,  and  P  values 
obtained  via  Wilcoxon  tests  for  nonsurviving  (NS)  versus  sur¬ 
viving  (S)  patient  groups  and  for  LSI  versus  non-LSI  (NLSI) 
patient  groups  are  shown  in  Tables  3  and  4,  respectively.  These 
tables  were  used  to  confirm  that  HRV  and  HRC  can  discrimi¬ 
nate  between  the  different  patient  groups,  respectively.  Heart 
rate  variability  (Poincare)  ratios  were  consistent  with  the  fact 
that  a  higher  ratio,  i.e.,  a  decreasing  long-term  variability  SD2, 
is  associated  with  increasing  risk  of  mortality  and  the  need  for 
an  LSI  (2,  14).  However,  only  maximum  HRV  ratios  for  NS 
patients  in  Table  3  were  statistically  larger  {P  <  0.05)  than  ratios 


Table  3.  Comparison  of  HRV  and  HRC  values  between  patient 
groups  for  mortality 


NS  patients  (n  =  1 7) 

S  patients  (n  =  91) 

P 

Measure 

Mean 

SD 

Mean 

SD 

Mean  HRV  ratio 

0.395 

0.198 

0.339 

0.208 

0.228 

Maximum  HRV  ratio 

0.927 

0.427 

0.713 

0.421 

0.063 

Minimum  HRV  ratio 

0.165 

0.146 

0.138 

0.072 

0.787 

Mean  HRC 

0.886 

0.343 

1.000 

0.310 

0.154 

Maximum  HRC 

1.549 

0.481 

1.591 

0.382 

0.524 

Minimum  HRC 

0.274 

0.275 

0.451 

0.348 

0.030 

SD,  standard  deviation. 


Table  4.  Comparison  of  HRV  and  HRC  values  between  patient 
groups  for  LSIs 


LSI  patients  NLSI  patients 

(n  =  82)  (n  =  26) 


Measure 

Mean 

SD 

Mean 

SD 

P 

Mean  HRV  ratio 

0.364 

0.230 

0.298 

0.091 

0.036 

Maximum  HRV  ratio 

0.811 

0.454 

0.545 

0.239 

<0.001 

Minimum  HRV  ratio 

0.133 

0.094 

0.172 

0.054 

0.011 

Mean  HRC 

0.908 

0.290 

1.217 

0.287 

<0.001 

Maximum  HRC 

1.535 

0.391 

1.743 

0.380 

<0.001 

Minimum  HRC 

0.321 

0.281 

0.744 

0.324 

<0.001 

for  S  patients.  Heart  rate  complexity  (SampEn)  values  for  pa¬ 
tients  who  received  at  least  one  LSI  were  consistent  with  the 
fact  that  this  group  often  has  lower  HRC  than  NLSI  patients 
(2,  14).  Only  minimum  HRC  values  for  NS  patients  in  Table  3 
were  statistically  smaller  than  the  values  for  S  patients. 

For  the  multivariate  logistic  regression  models  in  Table  5, 
results  showed  that  increasing  mean  HR  was  associated  with 
an  increased  risk  for  mortality.  Age,  height,  race,  and  weight 
were  removed  from  the  final  models  via  backward  elimination 
because  they  were  not  significantly  associated  with  mortality. 
In  the  model  for  vital  signs  alone  (Table  5),  the  odds  ratio  was 
1.04  (95%  confidence  interval  [Cl],  1.01-1.09;  P  =  0.007)  for 
mean  HR  (per  beats  per  minute  increase).  Importantly,  inclu¬ 
sion  of  HRC  in  the  multivariate  logistic  regression  analyses 
showed  that  decreasing  minimum  HRC  was  also  associated 
with  an  increased  risk  for  mortality,  increasing  odds  by  100%. 
Here,  odds  ratios  were  1.05  (95%  Cl,  1.01-1.09;  P  =  0.007)  for 
mean  HR  (per  beats  per  minute  increase)  and  0.01  (95%  Cl, 
0.00-0.58;  P  =  0.02)  for  minimum  HRC  (per  unit  increase). 

On  the  other  hand,  for  the  regression  models  in  Table  6, 
no  traditional  vital  signs  were  significantly  associated  with  an 
increased  risk  for  LSIs.  Again,  age,  height,  race,  and  weight 
were  removed  from  the  final  models  via  backward  elimination 
because  they  were  not  significantly  associated  with  LSIs.  In 
the  model  for  vital  signs  alone  (Table  6),  the  odds  ratio  was 
1.01  (95%  Cl,  0.99-1.04;  P  =  0.21)  for  mean  HR  (per  beats  per 
minute  increase).  Inclusion  of  HRC  in  the  regression  analyses 
again  showed  that  decreasing  minimum  HRC  was  associated 
with  an  increased  risk  for  LSIs  and  increased  odds  by  nearly 
100%.  The  odds  ratio  was  0.02  (95%  Cl,  0.00-0.14;  P  < 
0.0001)  for  minimum  HRC  (per  unit  increase). 

When  total  CCS  scores  were  included  in  the  models,  no 
traditional  vital  signs  were  significantly  associated  with  an 

Table  5.  Logistic  regression  models  with  various  risk  factors 
(excluding  GCS)  for  mortality 


Variable  Odds  ratio  for  mortality  (95%  Cl)*  P 


Field  HR 

1.04  (1.01-1.09) 

0.007 

With  HRC 

Field  HR 

1.05  (1.01-1.09) 

0.007 

Minimum  HRC 

0.01  (0.00-0.58) 

0.02 

*Odds  ratios  for  measurements  reflect  per-unit  increase. 
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Table  6.  Logistic  regression  models  with  various  risk  factors 
(excluding  GCS)  for  LSIs 

Variable 

Odds  ratio  for  LSIs  (95%  Cl)* 

P 

Field  HR 

1.01  (0.99-1.04) 

0.21 

With  HRC 

Field  HR 

1 .00  (0.98-1 .03) 

0.86 

Minimum  HRC 

0.02  (0.00-0.14) 

<0.0001 

*Odds  ratios  for  measurements  reflect  per-unit  increase 


increased  risk  for  mortality  or  LSIs.  However,  inclusion  of 
HRC  in  the  multivariate  logistic  regression  analyses  showed 
that  decreasing  minimum  HRC  was  also  associated  with  an 
increased  risk  for  mortality  or  LSIs  and  could  be  used  for  risk 
stratification.  Odds  ratios  for  regression  models  with  various 
risk  factors  (including  GCS  scores)  for  mortality  and  LSIs  are 
shown  in  Tables  7  and  8,  respectively. 

Receiver-operating  characteristic  curves  (Figs.  1  and  2)  de¬ 
monstrated  better  prediction  for  mortality  and  LSIs  using  HR 
and  HRC  (area  under  the  curve  [AUC];  AUCs,  0.86  and  0.86, 
respectively)  than  using  HR  alone  (AUCs,  0.79  and  0.57,  re¬ 
spectively).  Likewise,  ROC  curves  (Figs.  3  and  4)  demon¬ 
strated  better  prediction  for  mortality  and  LSIs  using  total  GCS 
score  and  HRC  (AUCs,  0.82  and  0.97,  respectively)  than  using 
total  GCS  score  alone  (AUCs,  0.81  and  0.91,  respectively). 

DISCUSSION 

This  study  investigated  the  effectiveness  of  using  traditional 
and  new  vital  signs  for  predicting  mortality  and  the  need  for 
LSIs  in  prehospital  trauma  patients  through  different  compar¬ 
isons  of  statistical  models,  ROC  curves,  and  AUC  results.  It 
was  the  first  to  use  a  multivariate  logistic  regression  model  to 
determine  whether  standard  vital  signs  along  with  HRV  and 
HRC  were  better  in  predicting  mortality  than  standard  vital 
signs  alone.  In  the  statistical  analyses  previously  described, 
HR  generally  increased  the  odds  of  death  or  an  LSI  by  no  more 
than  5%.  When  HRC  values  were  included  in  the  analyses, 
they  increased  the  odds  of  death  and  the  need  to  perform  LSIs 
by  approximately  100%.  When  GCS  scores  were  included  in 
the  analyses,  HRC  continued  to  increase  the  odds  of  death  and 
the  need  for  an  LSI  by  100%.  Following  recent  work  described 
in  Liu  et  al.  (2),  this  study  confirmed  that  models  using  both 
traditional  and  new  vital  signs  were  superior  over  models  using 
only  traditional  vital  signs  for  predicting  the  need  for  LSIs  in 
trauma  patients.  Furthermore,  this  study  confirmed  that  models 

Table  7.  Logistic  regression  models  with  various  risk  factors 

(including  GCS)  for  mortality 


Variable  Odds  Ratio  for  Mortality  (95%  Cl)*  P 


Total  GCS  score 

0.77  (0.62-0.90) 

0.001 

With  HRC 

Total  GCS  score 

0.79  (0.63-0.93) 

0.003 

Minimum  HRC 

0.20  (0.01-1.98) 

0.18 

*Odds  ratios  for  measurements  reflect  per-unit  increase. 


Table  8.  Logistic  regression  models  with  various  risk  factors 
(including  GCS)  for  LSIs 


Variable 

Odds  ratio  for  LSIs  (95%  Cl)* 

P 

Total  GCS  score 

0.41  (0.13-0.66) 

<0.0001 

With  HRC 

Total  GCS  score 

0.29  (0.04-0.63) 

<0.0001 

Minimum  HRC 

0.01  (0.00-0.10) 

<0.0001 

*Odds  ratios  for  measurements  reflect  per-unit  increase. 


using  both  traditional  and  new  vital  signs  were  superior  over 
models  using  only  traditional  vital  signs  for  predicting  mor¬ 
tality  in  trauma  patients.  The  major  implication  of  this  study 
was  that  traditional  and  new  vital  signs  should  be  used  simulta¬ 
neously  to  improve  prediction  of  mortality  and  the  need  for  LSIs 
in  prehospital  trauma  patients  during  all  echelons  of  trauma 
care.  Again,  this  study  demonstrated  that  multivariate  logistic 
regression  models  incorporating  HRV  and  HRC  could  increase 
the  mortality  and  LSI  prediction  accuracy  for  trauma  cohorts. 

Although  values  and  significances  of  odds  ratios  were  smaller 
than  those  described  in  Liu  et  al.  (2)  (most  likely  owing  to  the 
number  of  LSI  patients  versus  NLSI  patients),  this  study  pro¬ 
vided  a  different  cohort  to  validate  the  hypothesis  and  results 
presented  in  Liu  et  al.  (2).  It  is  important  to  note  that  this  study 
and  previous  work  (2)  involved  2  different  databases  containing 
disparate  sets  of  patients,  yet  both  studies  were  able  to  demon¬ 
strate  the  power  of  vital  signs  and  HRC  to  predict  the  need  for 
LSIs  in  prehospital  trauma  patients.  A  novelty  of  this  study  was 
the  exploration  of  traditional  and  new  vital  signs  for  predicating 
mortality,  showing  that  a  multivariate  logistic  regression  model 
combining  HR  and  HRC  was  superior  in  performance  over  a 
model  using  HR.  Unlike  the  methodology  described  in  Liu  et  al. 


1 -Specificity 
Survived 

Fig.  1 .  Receiver-operating  characteristic  curves  for  models  predicting 
mortality  (excluding  GCS  scores).  Receiver-operating  characteristic  curves 
were  obtained  to  examine  the  discriminating  power  of  multivariate  logistic  re¬ 
gression  models  (excluding  GCS  scores)  for  the  outcome  of  mortality  in  1 08 
subjects.  The  curves  demonstrated  better  mortality  prediction  for  models  using 
both  vital  signs  and  HRC  (AUC,  0.86)  than  for  models  using  only  vital  signs 
(AUC,  0.79). 
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Fig.  2.  Receiver-operating  characteristic  curves  for  modeis  predicting 
interventions  (exciuding  GCS  scores).  Receiver-operating  characteristic 
curves  were  obtained  to  examine  the  discriminating  power  of  multivariate  logistic 
regression  models  (excluding  GCS  scores)  forthe  outcome  of  at  least  one  LSI  in 
1 08  subjects.  The  curves  demonstrated  better  LSI  prediction  for  models  using 
both  vital  signs  and  FIRC  (AUC,  0.86)  than  for  models  using  only  vital  signs 
(AUC,  0.57). 
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Fig.  4.  Receiver-operating  characteristic  curves  for  modeis  predicting 
interventions  (including  GCS  scores).  Receiver-operating  characteristic  curves 
were  obtained  to  examine  the  discriminating  power  of  multivariate  logistic  regres¬ 
sion  models  (Including  GCS  scores)  for  the  outcome  of  at  least  one  LSI  in  1 08 
subjects.  The  curves  demonstrated  better  LSI  prediction  for  models  using  vital 
signs,  GCS  scores,  and  HRC  (AUC,  0.97)  than  for  models  using  only  GCS  scores 
(AUC,  0.91). 


(2),  HRV  was  calculated  via  the  Poincare  ratio  (see  Heart-Rate 
Variability  and  Complexity).  Moreover,  HRV  and  HRC  were 
retrospectively  derived  from  the  patients’  ECG  waveforms  rather 
than  extracted  from  the  vital  signs  monitor.  Owing  to  previous 
work  (2)  and  the  power  of  HRC  for  group  discrimination,  Tables  1 
through  4  only  showed  results  involving  HRC  rather  than  HRV 
(which  yielder  lower  odds  ratios  than  the  former).  Neverthe¬ 
less,  results  involving  HRV  were  similar  to  the  tabulated 


Fig.  3.  Receiver-operating  characteristic  curves  for  modeis  predicting 
mortality  (including  GCS  scores).  Receiver-operating  characteristic  curves 
were  obtained  to  examine  the  discriminating  power  of  multivariate  logistic  re¬ 
gression  models  (including  GCS  scores)  for  the  outcome  of  mortality  in  108 
subjects.  The  curves  demonstrated  better  mortality  prediction  for  models  using 
vital  signs,  GCS  scores,  and  HRC  (AUC,  0.82)  than  for  models  using  only  GCS 
scores  (AUC,  0.81). 


results.  For  example,  a  multivariate  logistic  regression  model 
combining  HR  values  and  maximum  Poincare  ratios  was  su¬ 
perior  in  performance  over  a  model  using  HR  alone,  yielding 
odds  ratios  of  1.06  (95%  Cl,  1.02-1.12;  P  =  0.002)  for  mean 
HR  (per  beats  per  minute  increase)  and  9.43  (95%  Cl, 
1.68-76.7;  P  =  0.011)  for  maximum  HRV  ratio  (per  unit  in¬ 
crease).  Inclusion  of  the  HRV  ratio  in  the  regression  analyses 
showed  that  increasing  the  maximum  Poincare  ratio  was  as¬ 
sociated  with  an  increased  risk  for  mortality  and  increased 
odds  by  nearly  10%.  Similarly,  ROC  curves  demonstrated 
better  prediction  for  mortality  and  LSIs  using  HR  and  HRV 
(AUCs,  0.86  and  0.73,  respectively)  than  using  HR  alone 
(AUCs,  0.79  and  0.57,  respectively). 

In  addition  to  the  aforementioned  notes,  it  is  important  to 
point  out  that  models  using  and  not  using  total  CCS  scores 
were  considered  separately  in  this  study  because  CCS  scores 
are  not  always  available  or  practical  to  obtain  in  prehospital 
and  battlefield  environments.  Like  previous  work  and  results 
described  in  Liu  et  al.  (2),  this  study  showed  that  models  not 
using  CCS  scores  were  still  able  to  predict  mortality  and  the 
need  for  LSIs  with  greater  than  80%  accuracy.  By  integrating 
other  mechanisms  for  scoring  injuries  in  a  rapidly  changing 
environment  (8),  it  is  possible  to  preclude  use  of  CCS  scores 
while  increasing  prediction  accuracy. 

Because  HRV  and  HRC  can  be  calculated  rapidly  and  non- 
invasively  (16,  20,  21),  they  may  be  leveraged  with  routine  vital 
signs  to  help  change  the  paradigm  in  prehospital  care  for  the 
severely  injured.  With  today’s  advances  in  monitoring  tech¬ 
nology,  these  parameters  could  be  displayed  simultaneously  on 
one  screen  and  used  at  the  injury  scene  for  risk  stratification 
of  mass  casualties  requiring  air  transport,  such  as  casualties 
occurring  on  the  battlefield  or  a  major  highway.  In  addition, 
standard  and  new  vital  signs  could  be  used  together  for  rapid 
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triage  of  prehospital  patients  (e.g.,  which  facility  or  what  level/ 
role  of  care)  or  during  transport  to  assist  with  timely  and 
accurate  administration  of  LSls.  Thus,  these  data  could  play  a 
role  in  facilitating  prehospital  care  and  patient  stabilization  until 
further  care  at  a  facility  is  available. 

Future  studies  may  include  indicators  of  numeric  and  wave¬ 
form  data  quality  to  provide  a  more  comprehensive  model  for 
predictions  of  outcomes  in  trauma  patients  as  well  as  development 
of  a  novel  prehospital  score  using  standard  vital  signs,  HRV, 
and  HRC. 


Limitations 

This  study  had  limitations  that  mirrored  previous  work  (2). 
In  addition,  lengths  of  ECG  waveforms  were  short,  and  no  data 
were  excluded  regardless  of  noise  in  the  captured  waveforms. 
A  third  limitation  was  that  manual  verification  of  all  R  waves 
for  108  trauma  patient  ECGs  was  performed  by  2  individuals 
having  no  background  in  cardiology.  Eourth,  this  study  did  not 
consider  analyses  for  examining  the  discriminating  power  of 
the  models  for  the  outcome  of  LSI  needs  resulting  in  mortality 
(i.e.,  patients  who  required  LSIs  and  who  eventually  died  after 
trauma).  A  strategy  similar  to  this  study  could  be  applied  to 
perform  these  analyses  in  the  future.  A  fifth  limitation  of  this 
study  was  that  times  of  actual  LSIs  were  not  recorded  and 
stored  in  the  TV  database  (8).  Hence,  this  study  could  not  di¬ 
rectly  measure  how  new  vital  signs  changed  with  LSIs  in  in¬ 
dividual  patients.  Nevertheless,  this  limitation  did  not  compromise 
the  results  of  this  study,  since  analyses  involved  extremes  (mini¬ 
mum  HRC  value  and  maximum  HRV  ratio),  thereby  strengthening 
the  potential  usefulness  of  new  vital  signs  for  triage  and  risk 
stratification. 

Here,  it  is  worth  clarifying  the  role  of  helicopter  medical 
transport  in  possibly  influencing  this  study’s  analyses  as  well 
as  the  limitations  of  HRV  and  HRC  in  the  clinical  settings. 
Assessment  of  breath  sounds,  bowel  sounds,  blood  pressure  by 
auscultation,  and  even  verbal  communications  with  the  patient 
is  greatly  impaired  during  helicopter  flight  due  to  noise  (5,  22). 
Eurthermore,  although  state-of-the-art  blood  pressure  monitors 
are  in  widespread  use  in  air  medical  services,  they  may  not 
always  be  reliable  at  low  pressures  (5).  Since  this  study  did  not 
focus  on  blood  pressures,  aforementioned  results  were  not 
compromised.  Inherent  in  this  study  is  the  possibility  of  me¬ 
chanical  noise  and  human  error  corrupting  segments  in  the 
ECG  waveforms,  which  could  have  altered  HRV  and  HRC 
calculations  for  analysis.  In  addition,  several  groups  (23,  24) 
have  reported  that  the  application  of  HRV  in  clinical  settings  is 
limited  by  poor  interindividual  and  intraindividual  reproduc¬ 
ibility.  Nevertheless,  their  conclusions  may  not  be  applicable 
to  all  HRV  or  HRC  measures  and  trends  in  these  measures. 
The  results  of  previous  work  (15)  and  this  study  have  shown 
that  even  in  the  presence  of  noise,  HRV  and  HRC  may  still  be 
able  to  discriminate  between  patient  groups  and  monitor  pa¬ 
tients  over  time.  Eurthermore,  recent  work  has  shown  that 
HRC  trends  can  be  reproduced  across  subjects  (25).  However, 
like  any  vital  sign,  HRV  or  HRC  alone  may  not  provide  a  com¬ 
plete  picture  of  patient  status.  Therefore,  new  vital  signs  such  as 


HRV  and  HRC  should  be  used  in  conjunction  with  standard  vital 
signs  and  other  parameters  for  improved  utility  (2). 

In  summary,  this  study  showed  that  multivariate  logistic 
regression  models  using  vital  signs,  HRV,  and  HRC  to  predict 
mortality  and  the  need  to  perform  LSIs  were  superior  in  pre¬ 
dictive  performance  over  models  using  solely  traditional  vital 
signs.  Improvements  in  the  timely  use  and  diagnostic  accuracy 
of  transportable  vital  signs  monitors  will  require  utility  of  tra¬ 
ditional  and  new  vital  signs  from  the  trauma  patient  cohort. 
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